Non-Markovian Diffusive Unravellings of Entanglement 
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The fully quantized model of two qubits coupled to a common bath is solved using the quantum 
state diffusion (QSD) approach in the non-Markovian regime. We have established an explicit 
time-local non-Markovian QSD equation for the two-qubit dissipative model. Diffusive quantum 
trajectories are applied to the entanglement estimation of two-qubit systems in a non-Markovian 
regime. In another interesting example, we have also considered exact entanglement unravellings for 
a dephasing model. In both cases, non-Markovian features of entanglement evolution are revealed 
through quantum diffusive unravellings in the qubit state space. 
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Many important realizations in quantum information, 
such as quantum computing, quantum communication 
and quantum cryptography, rely on the control and gen- 
eration of entanglement [1]. However, the true question 
arises in how to measure the entanglement of a quantum 
system in order to effectively use that information in ap- 
plication. For a quantum open system described by a 
reduced density matrix, most definitions of entanglement 
pertain to a property of an ensemble, such as entangle- 
ment of formation [2], E{p), and concurrence Q, C{p), 
however for possible applications in quantum informa- 
tion processing, the preparation of and measurement on 
a mixed entangled state would be quite cumbersome. A 
more approachable method would be to take advantage of 
the statistical nature of the quantum system and average 
over many realizations of a single system in order to infer 
information about the entanglement of the ensemble. Re- 
cently, entanglement unravellings in the Markov regime 
have been proposed 041] • For a general non-Markovian 
quantum open system, such a pure state approach is par- 
ticularly useful for the numerical simulation of the track- 
ing of entanglement information, which is known to be a 
hard problem due to the lack of a computable entangle- 
ment measure and a viable non-Markovian master equa- 
tion 0-[il- 

In this letter, our research serves as a first example 
of the efficient estimation of entanglement evolution in 
non-Markovian regimes without using the system den- 
sity matrix. We derive the exact quantum state diffusion 
(QSD) equation to estimate the entanglement evolution 
of a two-qubit system coupled to a bosonic heat bath at 
zero temperature [3-113 • For a general multi-qubit sys- 
tem, employing quantum trajectories over density matri- 
ces becomes enormously advantageous in terms of com- 
putational exhaustion. As will be shown in this letter, 
the entanglement computed from trajectories generally 
provides useful information about the status of the ac- 
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tual entanglement described by the system density ma- 
trix. For some initial states, the trajectory entanglement 
gives an almost identical estimation of the system entan- 
glement. 

We present a fully quantized model of two uncoupled 
qubits with respective transition frequencies uja and ujb 
that are coupled to a common zero-temperature heat 
bath via the interaction Hamiltonian H-mt = a Ifl i-^A + 
KLB)a'^ + g\{Ll^ + KL^g)a\], where La and Lb are Lind- 
blad operators describing the interaction of the qubits A 
and B with the heat bath, respectively, and k (0 < k < 1) 
is the control parameter that describes the ratio of the 
qubits' coupling strengths. The formal linear QSD equa- 
tion [l3| describing the dynamics of the quantum state 
of the qubits, tpt = tpt{x*), is given by 
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where the system Hamiltonian is i?sys = + "^^f 

and the Lindblad operator is £ = La + i^L b ■ The Gaus- 
sian process satisfies M[xtxl] — J2\ ISApe"*'^^^*"*) = 
a{t, s), the bath correlation function, where M[-] denotes 
the statistical mean over the noise. The solution to the 
QSD equation ipt, recovers the reduced density ma- 
trix of the qubit system: pt = M[\tpt){ipt\]- Central to the 
application of the QSD equation is to replace the func- 
tional derivative with a time-local operator, termed as 

0(t, s, x*)'ipt with initial 



the 0-operator, such that j^. 



condition 0{t — s,s,x*) — C. In principle, the existence 
of the 0-operator can be seen from the stochastic pro- 



pogator, liptix*)) = G{t,x*) |V'o) |2ll, but in practice it 
is difficult to find the explicit 0-operator. For the spe- 
cific two-qubit model presented in this letter, an exact 
equation for the 0-operator is derived upon satisfying 
the consistency condition (T3 |: 
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which ensures that ipt is a single-valued function, and 
thus establishes an exact solvable QSD equation. 
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We explore the non-Markovian regime by modeling the 
bath correlation function as an Ornstein-Uhlenbeck pro- 
cess such that a{t,s) = ^e"'''!*"^'. This continuous ran- 
dom process drifts toward a stationary long-term mean 
and is useful for viewing various memory effects via the 
parameter 7, which describes the rate at which noise that 
is progressing in time t becomes less and less correlated to 
its value at a particular past time s. As 7 grows very large 
and the correlation time Tc = ^ becomes very short, we 
will view the transition from non-Markovian to Marko- 
vian regimes and find that certain features are lost under 
the Markov approximation [l^-HEl- 

We investigate two types of qubit decoherence: dissi- 
pation of the quantum state and pure phase relaxation, 
both great challenges to maintaining robust entangle- 
ment. In both cases, we have derived the exact time-local 
0-operators, allowing us to efficiently solve the nonlinear 
QSD equation [T^, which robustly preserves the norm 
of the qubit state vector throughout every trajectory, a 
facet not exhibited by the linear equation [l^. The dy- 
namics of the normalized quantum state of the qubits, 
-ijjt = ji^Jfiji is presented as [l4 |: 

Jo 

-{{c^-{c^),)dit,s,i*))t]i't (3) 

where {A)^ — (iptlAl^t) is the quantum expectation value 
of operator A and x^ — x^ + dsa* {t, s) {C^) ^ is the 
shifted noise. Solving the stochastic differential equa- 
tion above for a particular realization of the Gaussian 
random noise reveals a single unraveling of the quantum 
system evolution, allowing one to calculate a single entan- 
glement trajectory by C(V't) = K^'tk^ ^ I'/'t )l- After 



a large number of realizations are produced, we take the 
mean over all concurrence trajectories, — M[C{tpt)], 
and obtain a value that is pertinent to the entanglement 
of the ensemble system. In this way, we can efficiently 
compute the approximate entanglement of a quantum 
open system without invoking the explicit form of the 
density matrix. The actual entanglement represented 
by the density matrix can be calculated through con- 
currence Q C{p) = max{0, v^Ai — y/X2 — \/X3 — V^} 
where {i = 1,2,3,4) are the eigenvalues of the ma- 
trix g = Pif^y ^ <^y)p*{'^y ^ <^y ) i^ desccndiug order. 
Upon direct comparison, it is clear that C,/, cannot be 
the same as true entanglement C (p) due to the concavity 
of the concurrence. However, can be used as an up- 
perbound of the actual entanglement, such that if « 
then C{p) « 0. In fact, as shown below, provides al- 
most perfect estimation of the actual entanglement for 
some initial states. Above all, the calculation of is 
much simpler than that of C(p), especially for a large 
number of qubits. This pronounces to be a good in- 
dicator for the actual behavior of the entanglement and 
will be explored in the upcoming models. 

A dissipative interaction, which causes the quantum 
state to lose energy, is denoted by the Lindblad op- 
erators La — (J^ and Lb = <J^ such that £ = 
+ Ka^. For the very first time, we have estab- 
lished an exact QSD equation for the dissipative inter- 
action model. By the consistency condition of Eq. ©, 
we find the exact 0{t,x*) = dsa{t, s)0{t, s,x*) to be 
6{t, X*) = A{t)a^ + B{t)cr^ + F{t)afcr^ + G{t)afa^ + 

i ^ J* ds' P{t, s')x*^ a^a^ , which is true for an arbitrary 

bath correlation function. By imposing the Ornstein- 
Uhlenbeck bath correlation function, we derive a set of 
differential equations for the coefficients of the O opera- 
tor: 



dtA{t) = -jA{t) + ^+ iojAA{t) + + 2KF{t)G(t) + G{tf - ^iQ{t), 

dtB{t) = --fB{t) + ^+iUBB{t) + KB{t)^ + 2F{t)G{t) + KF{tf-^iQ{t), 

dtF{t) = --fF{t) + iujBF{t) + F{t)[A{t) + G{t)] + B{t)[G{t) - A{t)] + 2KB{t)F{t) - \iQ{t), 

dtG{t) = --fG{t) + iujAG{t) + KF{t)[A{t) + G{t)] + KB{t)[G{t) - A{t)] + 2A{t)G{t) - qiQ{t), 

dtQ{t) = ~2-fQ{t) + i{ujA+ujB)Q{t) + 2[A{t) + KB{t)]Q{t)~i-i[F{t) + KG{t)\ (4) 



together with the explicit solution P(t,s') = —2i[F{t) + 
KG{t)] exp{/^* ds[-7 + iujA + iws + 2A(s) + 2nB{s)] } and 
initial conditions A{Q) = B{0) = G(0) = F(0) = Q(0) = 
0. It should be noted that when n — I and loa ~ l^b, 
A{t) = B(t) and F{t) = G{t), representing a highly 
symmetrical setup where qubit A and qubit B are in- 
terchangeable. Throughout the rest of the paper, we will 



assume loa = ojb = i^- 

Knowledge of the exact equations for the 0-operator 
allows us to solve the nonlinear QSD Equation for var- 
ious unravelings of the time evolution of the qubits ini- 
tially in the maximally entangled Bell States, l^**) = 
-^{ItDzblll)} for qubits with correlated spins and 

= ^{|tl> ± lit)} for qubits with anti-correlated 
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spins. In the interest of direct comparison to the ex- 
act case, we have derived the O operator in the Post- 
Markov approximation to be OpM — ifoi't)+i^fiit))C — 
h{t){af + «:2af )£, where h{t) = J^ait,s)d.s, h{t) = 

alt^sY^— s)ds^ and f2{t) = a{t, s)a{s,u){t — 

s)duds [22| . 




FIG. 1. Dissipative Model: The exact is compared to 
ensemble calculations in the non-Markovian regime, C{p), and 
trajectory methods under the Post-Markov approximation, 
CpmW with K = 1 and 7 = 0.3 for (a) |V'o) = |*+> and (b) 
l^o) = 1'!'+) 
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FIG. 2. Dissipative Model: For |V'o) = |*+>, over 1000 
realizations is compared over (a) various values of k for fixed 
7 = 0.3 and (b) various values of 7 for fixed k = 1. 
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FIG. 3. Dissipative Model: For It/jq) = |$^), is compared 
over (a) various values of k for fixed 7 = 0.3 for long times 
and (b) various values of 7 for fixed k = 0.25. 

The entanglement derived from the exact mean trajec- 
tory method, C'^ is compared to that under the Post- 
Markov approximation, CpM{ip), as well as the entan- 
glement of the ensemble, C(p), in Fig. [T] For both ini- 
tial states, 1*+) in (a) and 1$+) in (b), the actual en- 
tanglement C {p) displays the repeated revival and death 
of entanglement known to the two-qubit model j26 28]. 
Such trends are exhibited by the exact mean entangle- 
ment trajectory C^, whereas in contrast, applying the 



Post-Markov approximation removes all revival features 
of the curve. This figure demonstrates the dependence of 
the theory on the initial qubit state. In both cases, 
acts as an upperbound for the exact entanglement C{p), 
however |$^) provides a much closer approximation than 
1^'+). Upon taking a closer look at for various val- 
ues of K and 7 in Fig. [5] and Fig. [3J many interesting 
attributes of this model are revealed. 

For initial state \^o) = |*+), Fig. [DJa) highlights the 
significant revival feature of the equal couplings case, 
AC = 1 , which shrinks as we decrease the coupling strength 
of qubit B. On the other hand, the asymmetrical cou- 
plings case, K 7^ I, causes the entanglement to decay 
much slower. In this Non-Markovian region of 7 = 0.3 we 
find the entanglement will remain non-zero for extended 
time. In Fig. [Ifb) we look at a fixed k — 1 and view the 
transition from non-Markovian to Markovian regimes as 
we let 7 become large. It is shown that the revival peak 
grows as we tend toward the non-Markovian conditions 
and eventually oscillates very close to an entanglement 
state that will not decay, allowing one to maintain a 
highly entangled state over a long period of time when 
large memory effects are present. 

In Fig. 131a) we look at long time entanglment evolu- 
tion from initial state |$^) for various values of k and 
fixed 7 = 0.3 where many revival peaks are witnessed. 
Similar to the previous case of initial state IvP^), tall 
revival peaks are displayed for symmetrical couplings, 
however they come at the expense of a faster disentangle- 
ment. Once again, the k — case reveals a much slower 
entanglement decay and remains non-zero even for long 
times. Comparing the effects of memory on the entangle- 
ment dynamics, Fig.[31Jb) again demonstrates that a very 
large correlation time allows the quantum state to remain 
highly entangled for extended times. An interesting dif- 
ference for this initial state is that even for fairly large 
7 = 5, the rebirth of entanglement is still a dominant 
feature. 

As another important case, we consider a dephasing 
type of interaction, which provides an example of pure 
decoherence without dissipation. Described by the two 
Lindblad operators La = erf and Lb = erf , the consis- 
tency conditions of Eq. ([2]) result in the noise-free exact 
0-operator 0{t,s) — C — + Kaf . The mean entan- 
glement trajectories for the dephasing model are plotted 
in Figs. dH) and 

In Fig. lU^a), the dotted curve, k — 0, represents the 
scenario of qubit A interacting with the heat bath while 
qubit B is a free particle. As is expected of the single- 
qubit dephasing channel p^ . the entanglement of the 
qubits asymptotically decays to zero. Moreover, as we 
introduce the interaction of qubit B to the environment 
through K ^ 0, the disentanglement rate between the 
qubits only increases and causes a faster death of en- 
tanglement. In Fig. m^b), the very non-Markovian case, 
7 — 0.01, where the memory of the system extends much 
further into the past, reveals the preservation of entangle- 
ment for a considerable length of time before beginning 
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FIG. 4. Dephasing Model: For \'4>o) = 1*"^), over 1000 
realizations is compared over (a) various values of k for fixed 
7 = 1 and (b) various values of 7 for fixed k — 1. 
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FIG. 5. Dephasing Model: For |?/)o) = l^"^), is compared 
over (a) various values of k for fixed 7 = 1 and (b) various 
values of 7 for fixed k = 0.25. 



to decay. In the limit as 7 approaches one would ex- 
pect entanglement to be sustained at the maximum value 
eternally. As we shorten the memory of the system in the 
Markov regime, the entanglement curves reveal a steeper 



and steeper descent toward zero entanglement. 

The same analysis was applied to qubits with initial 
state 1$^) and shown in Fig. [S] In Fig.jSja) we immedi- 
ately see that, when the coupling constants of the qubits 
to the heat bath are equal, k — 1, then the initially en- 
tangled state is protected due to the symmetry between 
the two qubits. When the qubits are not coupled to the 
modes of the heat bath in exactly the same way, k 7^ 1, 
the entanglement will eventually decay to zero. Simi- 
lar to the dissipative model, Fig. [SJb) displays the pro- 
longed entanglement of the qubits in the non-Markovian 
case, 7 = 0.01, and the faster disentanglement rate of the 
Markov approximation, 7 = 10. 

In conclusion, we have shown that the dynamical en- 
tanglement of a non Markovian open system can be effi- 
ciently estimated by employing exact quantum diffusive 
trajectories. In particular, we have shown that entangle- 
ment dynamics of the system are very sensitive to which 
initial state the qubits evolved from, how the qubits are 
coupled to the heat bath, and the correlation time be- 
tween the noise variables. We emphasized that under 
the Markov approximation, the entanglement for both 
sets of Bell states was characterized by fast disentan- 
glement and suppressed revival features, whereas in the 
non-Markovian regime, large revivals were witnessed and 
for an extensively long correlation time, the qubits re- 
mained nearly maximally entangled for long times. This 
again emphasizes the importance of memory effects on 
the dynamics of a quantum system. 
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HROOll-09-1-0008 and the NSF PHY-0925174. 
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